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Wing flapping with minima energy 
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The interesting article by Archer* Sapuppo, and Betteridge^ leads 
to the problen of deternining the optlnun lift distribution along the wing 
during flapping notion. Adopting the assumptions made in Ref. 1* the 
problem may be reduced to one of minimizing the induced drag for a 
specified* periodically varying, bending moment at the wing root. The 
wing-root bending moment should be of primary importance to a bird or* 
for that matter* to a mechanical flapping machine, since it is at this 
point that most of the flapping energy is introduced. 

To solve this problem we adapt a technique employed by Prandtl and 
Hunk. Assume at first that no net lift is specified* that is* the wing 
is used solely for propulsion. We seek to determine the load distribu- 
tion £.(y) that minimizes the vortex drag with a given bending moment 
M. Assume that i(y) is given. Then small variations in the shape, 

6i(y) » will produce no first-order change in or M. Since two condi- 
tions are involved* it is sufficient to consider two discrete elements 
of variation (see Fig. 1). The added drag will have three components. 

(1) The drag of the added distribution , 6&2 in the downwash 
w^ of the original load distribution Jt(y). 
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(2) The drag of the original distribution in the downwash of the 
added distribution. 

(3) The drag of the added distribution alone. 

According to Hunk's reciprocal theorem (the "mutual drag" theorem), 
items (1) and (2) are equal; hence we need consider only item (1). 

Item (3) is proportional to the square of the added lift and hence does 
not contribute to the first-order variation. Ve then have the equations 

$*171 + $* 2 Y 2 “ 0 (for “ 0) (1) 

62^ + 5£ 2 w l2 - 0 (for GD ± - 0) (2) 

These equations will be satisfied for all positions of 5 £i, 622 ** we 
make 

w i " “i y * ( 3 ) 

where is a constant. Thus, we see that for minimum drag in flapping 
motion the induced downwash should vary linearly along the span; that 
is, the vortex trace of the wing should move as a rigid surface hinged 
at the wing root. 

In Hunk's original problem, the total lift and the wing span were 
given. There it was found that for minimum drag the induced downwash 
should be constant along the span or, in other words, the vortex trac^. 
of the wing should move downward as a rigid surface. This holds as 
well for wings having more complex shapes in front view, such as gull 
wings or wings with "winglets." Betz even extended this idea to 
propellers, showing that for minimum energy loss the vortex trace of a 
propeller should move backward as a rigid helix; the potential flow 
solution for this motion was found by Goldstein. That such ideas 
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should persist in the case of flapping notion is therefore not entirely 
unexpected. 

The extension of eqn. (3) to more complex shapes, such as gull 
wings, is easily accomplished by introducing the appropriate cosines 
into eqns. (1) and (2). In such cases it is found that minimum energy 
loss or drag occurs when the wake moves so as to satisfy the boundary 
condition of an impermeable surface having the shape of the wing trr :e 
and executing a similar motion. 

In eqn. (2) we assumed that the vortex drag could be calculated by 
considering an induced downwash w^ at the wing. This assumption limits 
our calculation to wings of high aspect ratio and to cases in which the 
wave length of the flapping motion is long compared with the wind chord. 
A still simpler situation arises if we consider very slow flapping such 
that the wavelength is large compared with the wing semispan (as assumed 
in Ref . 1) . In this case , the motion of the wake is approximately two- 
dimensional and conventional induced-drag theory can be used. 

Adopting the latter assumption, we have to find the spanwise load 
distribution that corresponds to a linear variation of the induced 
downwash w^/V. Using the standard inversion formula of airfoil theory 
we obtain (see Ref. 2) 

r ■ I V 2 + 4 co,h ‘‘ ■jTf) (4 > 

for the ideal distribution of circulation r during slow flapping (see 
Fig. 2). Here s is the wing semispan; it will be seen from eqn. (3) 
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that the quantity b^i/V ia the downwaah angle w^/V at the wing tip. 

If AC^ ia the additional lift coefficient during flapping we obtain 

«i» 3 

-r*jar“L • < 5 > 

If we let (d be the angular velocity of the flapping notion then' toy 
will be the vertical velocity of the notion at y. The instantaneous 
thrust will be 


T 


r 

l JoT(a - Mj] 


)y dy 


( 6 ) * 


Averaging over a cycle and forning the coefficient 




S | V 2 


we find 


^ (<0 - tOj)s WjS 


where to and to^ are now maximum values. The efficiency is 


to. 

_ _ JV_ . _i 

n 2Mb) ui 


(7) 


( 8 ) 


(9) 


where M is the bending moment. The loss of efficiency is thus sinply 
the ratio of the angular velocity induced by the wake to the angular 
velocity of the wing. At its maximum the efficiency is constant along 
the span. Introducing eqn. (9) into (8) we find 



i nU - n)(f ) 2 


( 10 ) 
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which is a simple equation connecting the thrust, the efficiency, and the 
tip helix angle. Figure 3 shows the relation between thrust and effi- 
ciency. Figure 4 shows the efficiency given by the optimum loading 
(eqn. (4)) compared with that determined in Ref. 3 for an elliptic load- 
ing. It will be noted th8t for a given value of thrust there are two 
possible values of the efficiency, the lower value evidently correspond- 
ing to an excessively high lift coefficient. 

Expressing in terms of the additional lift (eqn. (5)) we find 


n 


1 


_ 3 ^ 

4A tos/V 


(ID 


which shows the influence of aspect ratio. Integration over a periodic 
cycle shows that the efficiency is a maximum when the incremental lift 
AC^ is in phase with the upward or downward velocity of the wing tip. 
Hence, AC^ and id in eqn. (10) may be maximum values during the cycle. 
However, for best efficiency the lift increment should be small and the 
vertical velocity of the wing tip large-leading to rather large angles 
of twist during flapping, as found in Ref. 1. 

The exact angle of twist required to obtain the optimum loading 
will depend on planform of the wing. To compute this we make use of the 
formulas 

r - | c t v c , (12) 


where Cj is now the local lift coefficient and c the local wing chord. 
In addition we have 


C 


l 





+ Aa 



(13) 


where Aa 


g 


is the geometric twist angle introduced during the flapping. 
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Figure 5 shows curves of a /(ws/V) for elliptic wings of * 6 and 10. 
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It appears that the optimum twist for the elliptic wing is almost linear, 
conforming to the flapping motion. However, a small additional positive 
angle of attack appears to be needed at the wing root. This variation 
will, of course, be influenced by the planform shape of the wing. 

A complete analysis of the motion of a bird would involve vertical 
motions of the body in addition to the pure flapping motion considered 
here. Such vertical motions will not only Influence the optimum twist 
distribution, but will also contribute directly to the propulsion-possibly 
in a negative sense. 

Since the flapping wing is also normally used for support, it is 
necessary to consider the possible "interference drag" between the 
steady, supporting lift and the alternating lift due to flapping. 

According to Munk's reciprocal drag theorem we need to consider only 
one term of the interference drag. Taking the drag of the steady 
distribution in the downwash field of alternating distribution, it is 
easy to see that the interference is zero when averaged over a cycle, 
since the downwash of the added distribution is alternately positive 
and negative. 

It is of interest to compare the values of thrust given in Fig. 3 
with those required by a bird for steady flight. To maintain level 
flight we must have 


C„ - <L 


(14) 


The drag coefficient of a well streamlined large bird might be as low 
as 0.025. Assuming an aspect ratio of 8, 
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( 15 ) 



For a tip velocity ratio of 0.3, the propulsive efficiency is about 
87%. 

It has been suggested that the optimum load distribution during 
steady (nonflapping) flight would be the elliptic loading. However, as 
shown in Ref. 2, a 15% reduction of induced drag can be achieved for 
the same bending moment by adopting a more tapered form for the steady 
component of the load distribution. Although the variety seen in the 
wings of birds discourages generalization, this fact may explain why 
long-flying species, such as the gull and the Frigate bird, have more 
tapered wings. 

The foregoing°analysis is easily extended to the case of a propeller 
operating at a large advance angle. In this case, we are interested in 
the rolling moment or torque rather than the bending moment of the wing 
root. Equations (1) and (2) remain the same but the "downwash" given 
by eqn. (3) is now antisymmetric instead of symmetric and the absolute 
value sign on y is removed. The wake then revolves as a rigid helical 
strip, in agreement with the Betz-Goldstein result. 

The optimum spanwise loading in this case is simply 

T = 4 Ul s 2 (s]/i - ’ < 16 > 

giving the thrust coefficient, 

c t - 77n7> ' ? * " (1 ' n) ' (17) 
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where n is the efficiency. As in the case of flapping, the induced loss 
is the ratio of half the angular velocity of the far wake to the angular 
velocity of the wing or blade. 

Under the assumptions of slow flapping, or slow turning of the 
propeller, the normal slipstream loss considered in the actuator disc 
theory (Froude theory) will be negligible. The primary loss is the 
"rotational inflow" considered here. As is well known, this loss may 
be substantially overcome by using counterrotating propellers in tandem. 

The principle of the counterrotating propeller could be extended 
to wing flapping-a mode that might be termed "counterflapping." Two 
wings in a close tandem arrangement and moving in opposite phase would 
eliminate the induced losses we have calculated here. It is interesting 
that this mode of propulsion does not appear in nature-except possibly 
in the dragonfly (see Ref. 4). 
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FIGURE CAPTIONS 


Figure 1. Variational problem for flapping with minimum energy. 

Figure 2. Ideal distribution of the additional lift during slow flapping. 

Figure 3. Thrust and efficiency at different v&lues of the tip velocity 
ratio. 

Figure 4. Elliptic loading and optimum loading compared [v 1th] us/V 
- 0.3. 

Figure 5. Angle of twist during flapping elliptic planform: n ■ 80%. 
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Fig. 5 
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